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Chair:  Philip  L.  Marston 

~~)  A  scattering  algorithm  for  coated  spheres  was  developed  using  the  partial  wave 
series  of  Aden  and  Kerker  to  theoretically  determine  if  surface  films  on  air  bubbles  in 
water  have  a  significant  effect  on  optical  scattering  patterns.  The  air  bubbles  were 
modeled  as  spheres  (of  radius  a)  coated  by  a  film  of  uniform  thickness  h  and  complex 
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refractive  index  surrounded  by  water  of  refractive  index  n^  equal  to  4/3.  Size 

parameters  of  ka  (where  2ji/k  equals  the  optical  wavelength  in  water)  equal  to  100, 500, 
1000,  and  2500  (corresponding  to  radii  a  +  7.5  VMS)  to  189  p<n  for^t*^  ■  632.8  nm) 
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were  modeled  for  values  of  h  ranging  from  0  to  3  Jim.  Iiradiance  results  are  plotted  as  a 
function  of  scattering  angle  0  with  emphasis  cm  the  critical  scattering  region. 

Comparison  is  made  with  known  results  from  Mie  theory  for  noncoated  spheres  as  well  as 
a  physical-optics  approximation  developed  by  Marston  and  Kingsbury.  Some  modeling  is 
done  in  terms  of  geometric  ray  optics,  in  order  to  explain  observed  effects.  Finally, 
observations  regarding  the  use  of  noncoated  versus  coated  bubble  results  for  the  sizing 
and  detection  of  microbubbles  in  various  experimental  applications  is  discussed. 
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CHAPTER  I 


INTRODUCTION 


LL  Overview 

In  recent  years  a  great  deal  of  experimental  laboratory  research  and  computational 
modeling  has  been  done  regarding  the  behavior  of  light  scattered  by  air  bubbles  in  water.1'** 
In  these  investigations  the  bubbles  have  generally  been  modeled  as  a  non-absorbing 
uncoated  sphere  of  given  radius  a  and  relative  refractive  index  m  less  than  1.  Recent 
oceanographic  in  situ  experiments,  however,  have  suggested  that  bubbles  found  outside  the 
laboratory  may  possess  coatings  estimated  to  range  from  10  nm  to  1  pm  in  thickness.^' ^ 

If  this  is  indeed  the  case,  it  becomes  important  to  understand  and  model  the  effects  of  these 
surface  films  in  order  to  establish  the  range  of  validity  of  the  uncoated  laboratory  bubble 
results.  In  addition,  if  coated  bubble  results  were  found  to  vary  significantly,  the  modeling 
done  here  could  be  useful  in  determining  the  size  and  nature  of  the  coating  found  on  the 
bubble  which,  at  present,  is  still  unresolved. 

The  approach  taken  in  this  study  may  be  broken  down  as  follows;  first,  an  outline 
of  the  partial  wave  solution  used  for  the  case  of  a  coated  sphere  is  presented.  This 
mathematical  development  gives  some  insight  into  the  complexities  involved  in  the 
calculation  of  scattering  coefficients  and  also  serves  as  a  precursor  for  the  actual 
computational  results  which  follow  in  Chapter  2.  In  Chapter  2,  the  computational 
considerations  involved  in  writing  and  testing  the  coated  sphere  algorithm  are  discussed. 
The  results  for  the  range  of  parameters  tested  are  then  presented  along  with  a  discussion  of 
noted  effects.  A  comparison  is  also  made  with  Mie  results  for  a  noncoated  sphere  as  well 
as  a  physical-optics  approximation  developed  by  Marston  and  Kingsbury.2  Using  this 


The  exact  solution  to  the  problem  of  the  scattering  of  electromagnetic  waves  from 
concentric  spheres  was  first  presented  in  a  paper  by  Aden  and  Kerker  in  1951. 14  It  was 
based  on  a  development  presented  in  1908  by  Gustave  Mie  regarding  the  field  scattered  by 
any  spherically  symmetric  particle  of  arbitrary  radius  and  refractive  index  (commonly 
referred  to  as  Mie  theory).  Although  the  solution  to  the  problem  of  scattering  by  coated 
and  noncoated  spheres  has  been  available  for  many  years,  it  is  only  since  the  advent  of 
large  digital  computers  that  it  has  become  practical  to  use  in  scattering  The 

reason  for  this  will  become  apparent  in  the  following  development  for  an  electromagnetic 
wave  incident  on  a  coated  sphere. 

Consider  an  electromagnetic  plane  wave  which  is  incident  upon  a  coated  sphere  of 
radius  a  and  outer  radius  b  as  is  shown  in  Fig.  1.1.  The  relative  refractive  index  of  the 
core  is  mi  «  n/n,,  while  that  of  the  coat  is  =  r^/tiQ  where  no,  n,.  and  nj  are  the 

refractive  indices  of  the  outer  medium,  coating  and  inner  medium  respectively.  For  our 
particular  application  of  a  coated  air  bubble  in  water  no  *  nw,  the  refractive  index  of  water 
and  nj  *  n^  In  order  to  determine  the  fields  in  the  given  regions;  0sr£a,a<r<b 
(a  #  0)  and  r^b  it  is  necessary  to  solve  the  pertinent  vector  wave  equations  given  by 
Maxwell's  relations  using  spherical  geometry  and  the  appropriate  boundary  conditions. 
The  mathematics  required  to  do  this  is  straightforward  but  extremely  tedious  and  therefore 
steps  have  been  omitted  and  emphasis  placed  on  the  solutions.  The  reader  is  referred  to 


4 
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texts  by  van  de  Hulst15  and  Bohren  and  Huffman16  for  a  more  thorough  description.  The 
following  outline  was  taken  from  Chapter  8  of  the  Bohren  and  Huffman  text. 

Considering  again  the  situation  shown  in  Fig.  1.1,  the  solution  of  the  scalar  wave 
equations  in  spherical  polar  coordinates  yields: 

Verna  ■  coa  m*  P„m  (cos0)z„(kr),  (1) 


Vomn  =  sinm*  Pnm  (cos6)zn(kr),  (2) 

where  Pnm(cos0)  are  the  associated  Legendre  functions  and  is  any  of  the  four 
spherical  Bessel  functions  j„,  y„,  hn0),  or  h„(2)  depending  on  the  region  of  interest.  The 

subscripts  e  and  o  denote  even  and  odd.  Due  to  the  completeness  of  the  functions 
cos  m<p,  sin  mtp,  Pnm(cos0)  and  z„(kr),  any  function  that  satisfies  the  scalar  wave 

equation  in  spherical  polar  coordinates  may  be  expanded  as  an  infinite  series  in  the 
functions  (1)  and  (2).  The  vector  spherical  harmonics  generated  by  yemn  and  Vomn 
are 


^®emn  *  V  x  (r  Ventn)*  Momn  =  ^  x  (r  Vomn)» 

^cmn  *  (V  x  Menui)/k,  ^omn  ®  (V  x 

which  are  written  out  in  component  form  in  Bohren  and  Huffman. 17  Finally,  the  incident 
wave  is  taken  to  be  of  the  form  Ej  ■  Eq  e^r  *  ®*)0X  where  tx  indicates  the  x  unit 
vector  in  spherical  coordinates.  In  the  following  description  we  omit  the  exp(-icot)  time 
dependence  for  all  fields. 

With  this  information  it  is  now  possible  to  write  out  the  fields  for  the  various 
regions  of  the  coated  sphere.  The  electromagnetic  field  (Ej,  Hj)  in  the  region  from 
O^r^a  is  given  by 
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Ei  -  I  EJcBMolnO)-idnNeinO)l. 
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Hi  *  -kj/wm  I  E„[dnMeln(»  +  icnNolnO)lf 
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(3b) 


where  E„  *  inE0(2n+l)/n(n+l),  m  is  the  permeability  of  region  1,  ki  is  the  wave 
number  of  the  wave  with  frequency  to  in  that  region  and  c„,  dn  are  coefficients 
determined  by  the  boundary  conditions.  The  scattered  field  (Es,  Hs)  is  given  by 

E,  -  I  EnlianNemOl-bnMo^)],  (4a) 

n=«l 

H,  -  Wo*  I  +  (4b) 

n=l 


where  again  a„  and  bn  are  coefficients  which  will  be  determined  by  the  appropriate 

boundary  conditions.  The  (3)  superscript  indicates  that  the  outgoing  wave  solutions  in 
Eqs.  (1)  and  (2)  with  z„(kr)  =  ^(^(kjr)  are  to  be  used.  Due  to  the  requirement  of 

finiteness  at  the  origin,  the  radial  part  of  the  functions  in  Eqs.  (1)  and  (2)  which  generate 
the  vector  harmonics  in  die  expansion  of  (Et,  Hi)  is  constrained  to  be  jn  [indicated  by 
(1)  superscript].  However,  in  the  region  a  £  r  £  b  (a  #  0),  both  the  spherical  Bessel 
functions  j„  and  y„  are  finite;  as  a  consequence  the  field  in  region  2  (E2,  Hj)  is  written 
as 

OO 

E 2  *  n5j  Ejil^n  Motn(,)  *  »gn  Neln(1)  +  vn  Mqi^)  -  iwn  Nejn(2)],  (5a) 


H2  *  -  k2/o^2  I  En[gn  Meln(»  +  ifnN0in(')  +  w„  Mein<2>  +  ivnN0i„(2)],  (5b) 


where  the  vector  harmonics  Mcl„(2)  and  so  on,  are  generated  by  functions  of  the  form  of 
Eqs.  (1)  and  (2)  with  radial  dependence  y„(k2r).  The  boundary  conditions 


(E2  *  Et)  x  Sf  «  0,  (H2  -  Hi)  x$r  =  0,  r  =  a 

(Es  +  Ej - E2) x  £r  *  0,  (Hs  +  Hi-H2)xer  »  0,  r*b 
yield  eight  equations  in  the  coefficients  a,,,  b„,  c„,  d„,  fn,  gn,  v„,  wn: 

fn“lVn(®2*)  -  vnmlXn(m2x) '  Cr^VnC®^)  =  0,  (6a) 

wnmlX'n(m2x)  *  gnmiVn(m2x)  +  dnWn^l*)  -  0-  (6b) 

vnHlX’n(m2x)  *  fntilV'n^x)  +  CnH2Vn(miX)  *  0,  (6c) 

gnPl¥n(®2x)  -  wnp.1xn(m2x)  +  d^VnC®^)  -  0,  (6d) 

®2V'n(y)  *  anra2^n(y)  *  gnV’n^)  +  wnX'n(m2y)  =  0,  (6e) 

®2bn^n(y)  *  ®2Vn(y)  +  ^Vn(®2y)  *  v^m^)  =  0,  (6f) 

ti2Vn(y)*^24„(y)-gn^Vn(®2y)  +  wnllXn(®2y)  -  0,  (6g) 

bnM2^'n(y)  *  42¥*n(y)  +  fnPV’n<®2 y)  '  v#X’n(®2y)  *  0.  (6h) 


where  m|  and  are  the  refractive  indices  of  the  core  and  coating  relative  to  the 
surrounding  medium;  p,  P],  p2  a®  the  permeabilities  of  the  surrounding  medium,  core, 

and  coating;  and  x  ■  ka,  y  *  kb.  The  prime  indicates  differentiation  with  respect  to  x  or 
y  in  the  argument.  The  Riccati- Bessel  function  Xn(z)  *s  -*yn(z)-  Let  us  assume  for 
simplicity  that  p  =  Pi  =  p2  and  solve  the  set  of  above  equations  for  the  scattering 
coefficients  a„  and  b-:18 
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*"  "  §«i(y)[V'n(m2y)  -  AnX'n(“2y)]  ‘  “24'n(y)[¥n(m2y)  -  AnXn(m2y)]  ’ 

.  a  m2yn(y)[v‘n(m2y)  •  B^m#)}  •  ¥'n(y)[¥n(m2y)  -  BnXn(m2y)] 
m2§n(y)f¥n(m2y)  *  BnX'n^y)]  *  5’n(y)[¥n(m2y)  -  BnXn(m2y)] 


An 


m2¥n(m2X)¥,n(mix)  -  m1¥n(m2x)yn(m1x) 
m2Xn(m2x)¥n(ni l x)  -  miX'n^xtynOntx) 


(7c) 


B„ 


m2¥n(m  1  x)¥* n(m2x)  '  tPl¥‘n(m2x)¥*n(mlx) 
m2Xn(m2x)¥n(m  l x)  -  m1¥n(m1x)xn(m2x) 


(7d) 


Thus  we  find  that  the  determination  of  scattering  coefficients  a„  and  b„  becomes  an 
impressive  task.  However,  with  the  aid  of  a  computer  the  computation  can  be 
accomplished  as  will  be  shown  in  Chapter  2. 

For  future  reference  it  is  desirable  to  define  the  following  complex  scattering 
amplitudes  which  are  descriptive  of  the  scattered  fields:19 
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Sj(cos0)  =  X  2n^1  [anJtnfcosQ)  +  b„Tn(cos0)J, 
n=l  n(n+l) 


(8a) 


S2(cos0)  =  X  2n+*  [antjjfcosO)  +  bnJt„(cos0)]. 
n*l  n(n+l) 


(8b) 


The  angular  functions  Iq  and  k„  are  defined  by  the  relations 

Jtn(cos0)  =  — —  PnkcosO), 
sin0 


*n(cos0) 
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CHAPTER  2 


COATED  BUBBLE  SCATTERING  COMPUTATIONS 

2.L  Computational  Considerations 

A  computational  scheme  is  presented  in  this  section  for  the  calculations  of  the 
scattered  irradiances  Ij,  where  j  =  1  indicates  polarization  of  the  incident  electric  field 

vector  perpendicular  to  the  scattering  plane  and  j  =  2  indicates  the  parallel  case.  It  is 
based  on  similar  schemes  developed  by  Wiscombe1*2  and  Bohren  and  Huffman3  but 
differs  in  that  it  computes  the  coated  sphere  scattering  irradiance  as  a  function  of  scattering 
angle  rather  than  just  the  efficiency  factors.  Another  difference  is  that  it  also  uses  down 
recursion  for  all  logarithmic  derivative  calculations. 

Adopting  the  notation  of  Bohren  and  Huffman,  we  may  calculate  irradiances  (Ij) 
by  way  of  Eqs.  (8a)  and  (8b)  using  the  values  of  the  coefficients  a„  and  bn  determined 
in  Eqs.  (7a)  and  (7b).  The  logarithmic  derivatives  of  the  Ricatti-Bessel  functions  found 
in  Eqs.  (7c)  and  (7d)  were  calculated  by  the  down  recursion  technique  for  Mie  scattering 
developed  by  Wiscombe,  which  is  known  to  be  stable.  This  technique  consists  of 
specifying  a  priori  the  last  term  in  the  partial  wave  series  denoted  by  N  and  then  using 
recursion  relations  to  calculate  down  to  the  n  =  0  term.  By  doing  this  the  computational 
error  is  decreased  in  each  step.  The  problem  with  the  method,  however,  lies  in  selection 
of  the  N*  term.  Wiscombe,  through  extensive  testing,  has  determined  N  =  kb  +  4.05 
(kb) 1/3  +  2  (where  here  b  is  the  radius  of  the  outer  sphere  shown  in  Fig.  1.1)  is  the 
appropriate  relation  to  determine  N  and  hence  this  was  the  expression  used  in  our 
computations.  Another  necessary  consideration  was  the  preliminary  determination  of 
angle  step  size  N^g,  which  increases  in  proportion  to  ka.  In  order  to  resolve  the  fine 


structure  shown  in  Mie  results,  the  value  of  N^g  <  A/7  a  rad  (where  a  is  the  inner 

sphere  radius  shown  in  Fig.  1.1)  was  taken  from  previously  published  bubble  results.4 

The  partial  wave  scattering  coefficients  a„  and  bn  were  computed  using  the  IBM 

FORTRAN  VS  program  which  is  given  in  Appendix  A.  All  computations  were  performed 
on  an  IBM  3090-200  computer.  The  program  was  written  using  IBM  quadruple  precision 
(32  significant  digits)  as  this  insured  that  the  effects  of  the  coating  would  not  cause 
numerical  instabilities.  The  complex  arithmetic  was  done  using  IBM  FORTRAN  VS 
subroutines  again  using  quadruple  precision  which  allowed  both  the  real  and  imaginary 
components  in  all  calculations  to  have  32  significant  digit  precision. 

The  program  was  tested  in  a  wide  variety  of  ways.  Limiting  cases  were  taken  in 
which  the  coating  was  allowed  to  go  to  zero  thickness  which  gave  the  expected  noncoated 
bubble  Mie  results  previously  published.  3  1116  coating  was  also  given  die  ««n*»  index 
of  refraction  as  the  surrounding  medium  (water)  which  again  yielded  the  correct  Mie 
results.  Next,  results  were  checked  against  other  coated  sphere  calculations  done  by 
Cooper  et  al.®  and  B  runs  ting. ^  We  were  able  to  duplicate  Cooper’s  results  exactly  for  all 
cases  and  Brunsting's  for  scattering  angles  ranging  from  approximately  6  =  0°  to  90°. 
Scattering  coefficients  and  efficiency  factors  calculated  by  our  program  were  checked 
against  results  from  the  subroutine  written  by  Bohren  and  Huffman8  and  were  also  in 
perfect  agreement 

Due  to  the  use  of  quadruple  precision  and  the  need  to  calculate  three  logarithmic 
derivatives  by  down  recursion  the  program  is  rather  slow.  Typical  times  ranged  from 
approximately  10  c.p.u.  seconds  for  ka  «  100  with  a  1  )im  coating  to  roughly  150  c.p.u. 
seconds  for  ka  *  2500  with  a  coating  of  1  )ira.  Angle  step  size  was  typically  on  the  order 
of  (A/14  a)  rad  to  avoid  sampling  errors  for  any  size  coating  calculatic  •'ntrast  this 

the  physical-optics  approximation  program  (given  in  Appendix  A),  whic.  .e  to  predict 


the  coarse  structure,  executes  in  a  few  seconds  regardless  of  the  size  of  ka.  The  equations 
used  to  develop  this  program  are  listed  at  the  end  of  Section  II  of  Ref.  9.  Comparison  of 
coated  bubble  cases  will  be  made  with  the  physical-optics  approximation  for  forward 
scattering  (0  *  0°  to  90°)  to  determine  its  usefulness  in  coarse  structure  calculations  for 
coated  bubbles. 

2u2u _ Results 

A.  Effect  of  Increasing  Coating  Thickness 

A  typical  set  of  calculations  for  a  coated  air  bubble  in  water  with  ka  =  100 
(a  »  7.5  pm  for  »  474.6  nm)  over  the  0  ■  0°  to  180°  range  are  shown  in  Figs.  2. 1 
and  2.2.  The  parameters  used  were  the  following:  the  index  of  refraction  of  the  coat  tig 
was  real  (nonabsorbing  coat)  and  equal  to  1.5,  the  index  of  refraction  of  the  water  nw  wa: 
also  real  and  equal  to  4/3.  The  thickness  of  the  coatings  h  were  0.26  pm  and  1.14  pm 
respectively.  The  values  of  the  coating  were  determined  by  taking  the  ratio  of  the  inner 
sphere  radius  a  and  the  outer  sphere  b  (described  by  expression  AOB  in  the  COATSPHR 
program  in  Appendix  A).  The  coating  thickness  h  was  then  determined  by  the  relation 
h  =  b-  a- b(l  -  AOB).  The  range  of  bubble  sizes  can  be  illustrated  by  considering  HeNe 
laser  illumination  with  a  wavelength  in  air  equal  to  632.8  nm.  The  wavelength  in 
water  is  then  given  by  X*  ■  Xajl/nw  *  474.6  nm.  Using  X^  to  determine  a  for 
ka  «  100  gives  a  »  7.5  pm.  A  list  of  bubble  sizes  in  terms  of  h  and  the  actual  AOB 
values  used  in  the  calculations  for  X*  »  474.6  nm  are  given  in  Table  2.1.  Some  concern 
may  be  raised  about  using  the  approximate  value  of  nw  =  4/3;  however,  computational 
testing  over  a  range  of  nw  from  to  1.332  (water  at  30°C  for  X^  =  589.2  nm)  to  1.338 
(high  salinity)  showed  the  4/3  value  to  be  a  very  good  approximation.  Representative 
results  of  this  evaluation  are  given  in  Appendix  C.  Normalization  in  the  figures  was 
chosen  so  that  lj(0)  ■  1  represents  perfect  reflection  from  an  uncoated  sphere  of  radius  a 


Table  2.1 


Coating  Thickness  Values  for  X*  =  474.6  nm  (Xgj,  =  632.8  nm) 

This  *ahi#»  may  be  used  to  find  the  actual  value  of  AOB  used  in  the  computations  except  in 
certain  cases  where  the  value  of  AOB  is  noted  in  the  figure  caption. 


Figure  2.1.  Calculated  nonnalized  scattering  ixradiances  for  ka  »  100,  nw  »  4/3,  nc  a  1.5, 

0  »  0°  to  180°  for  the  log  of  the  perpendicular  polarization  Ij.  The  solid  curve  is  from  Mie 

theory.  The  thiixiashed  curve  is  the  calculated  coated  sphere  result  for  (a)  h  *  0.26  pm, 

(b)  h  «  1.14  pm.  All  values  of  h  are  for  HeNe  laser  light  in  water  with  ^  =  632.8  run 
The  inner  radius  a  •  7.5  pm. 
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Fig.  2.2  (a)  and  (b) 


according  to  geometric  optics.10  At  a  distance  R»ka2  from  the  center  of  the  bubble, 
the  actual  j-polarization  irradiance  is  incident  j-polarized  inadiance  multiplied  by 
IjKa/R)2/*]  where  Ij  -  4ISjl2  (ka)*2  with  j  again  denoting  the  previously  mentioned 

polarization  cases. 1 1 

Analysis  of  Mie  results  (indicated  by  solid  curve  in  Figs.  2.1  and  2.2)  has 
shown12  two  distinct  regions  of  interest  for  detecting  and  sizing  of  noncoated  bubbles  in 
water.  The  first  of  these  regions  known  as  the  Brewster  scattering  angle  region,  is  the 
region  around  the  Brewster  angle  determined  by  the  condition  0g  *  jc  -  2tan'1(m1)  (63  - 
106.3°  for  a  noncoated  air  bubble  in  water)  where  mi  *  n1j1/nw.  The  second  region  is 
that  near  the  critical  angle  given  by  0C  -  *  -  2sin*1(mi)  (0C  -  82.8°  for  noncoated  air 
bubbles  in  water).  The  critical  angle  region  irradiance  is  generally  used  for  size 
determination  while  the  ratio  of  the  two  regions  may  be  used  to  determine  if  die  scattered 
irradiance  is  indeed  from  a  bubble.  Furthermore,  data12  and  inspection  of  Mie  results  for 
uncoated  bubbles,  suggests  that  the  normalized  irradiance  for  the  j »  2  case  is  closer  to 
the  universal  value  of  Ij  ■  1/4  predicted  at  0C  by  the  physical-optics  approximation.^ 

For  this  reason  and  the  fact  that  parallel  polarization  has  a  simpler  fine  structure  the  j  =  2 
case  is  preferred  for  sizing  and  detecting  microbubbles  in  water.  Therefore  in  further 
illustrations  we  shall  use  the  j  -  2  case  while  the  supplementary  j  *  1  cases  are  given  in 
Appendix  B. 

Figures  2.1  and  2.2  show  that  even  with  the  coating,  the  coarse  and  fine  structure 
oscillations  described  by  Langley  and  Marston12  still  exist.  However,  with  increasing 
coating  thickness  some  structural  changes  become  apparent.  An  example  of  these  changes 
in  the  Brewster  angle  region  is  shown  in  Figs.  2.3(a)-(c)  and  2.4(a)-(c)  (thin-dashed 
curve).  Through  inspection  of  Figs.  2.3  and  2.4  one  finds,  in  general,  that  the  effect  of 
the  coating  is  to  raise  the  baseline  irradiance  near  0g  and  to  increase  the  magnitude  of  the 


*2-  ine  solid  curve  is  trom  Mie  theory  and  the  thin-dashed  curve  is 
sphere  result.  With  =  632.8  nm  the  coating  thickness  values 
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fing  structure  oscillations  (especially  prominent  in  Fig.  2.4(c)).  This  behavior  may  affect 
the  use  of  the  Brewster  angle  region  as  a  method  of  discriminating  bubbles  from  particles 
in  water  since  the  ratio  of  the  irradiances  between  it  and  the  critical  region  vary  non- 
uniquely  for  different  coating  values.  The  critical  angle  region  also  shows  this  behavior 
but  due  to  the  fact  that  iiradiance  values  are  generally  several  orders  of  magnitude  higher 
the  effects  axe  of  greater  interest  for  the  sizing  of  coated  bubbles  and  hence  further  analysis 
shall  be  concentrated  on  this  region. 

It  has  been  demonstrated  by  Mars  ton  and  Langley14  that  noncoated  bubbles  exhibit 
coarse  irradiance  oscillations  as  the  scattering  angle  6  decreases  below  a  critical  value  for 
total  reflection  6C  -  82.82°  in  water.  With  the  introduction  of  a  coating  under  the 
condition  n,.  >  nw,  we  find  that  this  coarse  structure  begins  to  shift  towards  this 

noncoated  critical  angle  with  the  degree  of  the  angle  shift  dependent  on  the  coating 
thickness  h  and  the  coating  index  of  refraction  n^  In  addition  to  this  coarse  structure 
angle  shift  the  magnitude  of  both  It  and  I2  can  generally  increase.  Examples  of  this 
angular  shift  and  magnitude  amplification  for  ka  =  100  with  thickness  values  ranging  from 
0.25  pm  to  3.0  pm  are  shown  in  Fig.  2.5(a)-(e).  (Note  that  in  these  figures  coarse 
structure  comparison  is  also  made  with  the  physical-optics  approximation.)  The  fact  that 
the  shifting  of  coarse  structure  is  small  near  6C  is  noteworthy  since  noncoated  results 

may  be  used  to  estimate  the  magnitude  of  the  scattering  from  coated  bubbles  in  many  cases 
of  interest  Note  also  that  the  physical-optics  approximation  does  very  well  for  coated 
bubbles  in  the  critical  angle  region  for  thickness  values  less  than  1  pm.  The  next 
structural  change  examined  is  the  change  in  magnitude  and  period  of  the  superposed  fine 
structure. 

Figure  2.6  illustrates  the  evolution  of  the  fine  structure  in  I2  from  30°  to  90°  as  a 
function  of  coating  thickness.  In  this  figure  the  effects  of  the  coating  can  be  seen  as  early 
as  h  *  0.08  pm  (Fig.  2.6(a)).  As  the  thickness  increases,  both  the  first  and  second  peaks 


Figure  2.5.  Calculated  normalized  irradiances  for  lea  =  100  in  the  critical  angle  scattering 
region  for  I2.  The  solid  curve  is  from  Mie  theory,  the  thin-dashed  curve  is  the  calculated 


coated  sphere  result,  and  the  thick-dashed  curve  is  from  the  Marston-Kingsbury  physical- 
optics  approximation.  The  coating  thickness  values  for  *  632.8  nm  are 

(a)  h  *  0.25  Jim,  (b)  h  *  0.5  Jim,  (c)  h  ■  1.0  Jim,  (d)  h  »  2.0  Jim,  and  (e)  h  *  3.0  Jim. 


pr-wj 


Fig  2.5  (e) 


Figure  2.6.  Evolution  of  the  fine  structure  as  a  function  of  coating  thickness  for  parallel 
polarized  irradiance  for  ka  -  100.  Coating  thickness  values  for  -  632.8  nm  are 

(a)  0  pm,  (b)  0.08  pm,  (c)  0.15  pm,  (d)  0.19  pm,  (e)  0.23  pm,  (f)  0.48  pm, 

(g)  0.66  pm,  (h)  0.93  pm,  and  (i)  1.13  pm.  The  corresponding  values  of  AOB  specified 
in  the  computation  were:  (a)  1.0,  (b)  0.990,  (c)  0.980,  (d)  0.975,  (e)  0.970,  (f)  0.940, 


from  the  left  begin  to  decay  into  finer  structure  while  the  large  number  of  smaller  peaks 
increase  in  amplitude  and  move  closer  together.  At  h  *  0.23  ^m  (Fig.  2.6(e))  the  fine 
structure  becomes  quite  prominent.  An  interesting  effect  was  also  noted  certain  values 
like  h  *  0.66  p.m  (Fig.  2.6(g)).  At  these  particular  values  the  fine  structure  amplitude 
is  diminished  as  is  shown.  Possible  reasons  for  this  will  be  discussed  in  the  next  section. 

B.  Geometric  Modeling 

Interpretation  of  the  results  found  in  Section  A  may  be  done  using  Fig.  2.7. 
Adopting  the  notation  of  Langley  we  define  rays  using  the  notation  (p,/),  where  p  is  the 
number  of  chords  within  the  bubble  and  /  equals  the  number  of  crossings  of  the  optic 
axis.  Note  that  a  subscript  b  will  be  used  to  indicate  rays  that  pertain  to  the  outer  sphere 
only.  Examining  the  two  rays  used  in  the  physical-optics  approximation,  namely  the  (0,0) 
and  (1,0)  rays,  and  the  (2,1)  ray  the  figure  shows  that  die  coating  shifts  the  scattered  ray 
towards  higher  values  of  0  as  well  as  increasing  the  distance  between  the  (0,0)  and  (2,1) 
ray.  The  first  mentioned  effect  partially  accounts  for  the  observed  coarse  structure  angle 
shift  while  the  second  partially  accounts  for  the  change  in  the  fine  structure. 

It  has  been  shown  that  interference  between  the  (0,0)  and  (2,1)  rays  is  partially 
responsible  for  the  fine  structure  oscillations.  The  quasi-period  A0fs  of  this  structure 

near  the  critical  angle  has  been  calculated  in  a  first  term  approximation  to  be  roughly  equal 
to  1VB2*  where  B2  is  the  distance  between  the  (0,0)  and  the  (2,1)  ray.  Thus,  by 
increasing  B2  with  the  addition  of  the  coating,  the  fine  structure  will  have  a  smaller 
penod  which  is  clearly  the  case.  The  damped  fine  structure  seen  in  Fig.  2.6(g)  is  most 
likely  due  to  effect  of  the  coating  on  reflection  and  transmission  coefficients.  The  purpose 
of  the  geometric  modeling  is  to  give  the  reader  some  physical  insight  in  the  interpretation 
of  the  computational  results  presented.  It  is  by  no  means  meant  to  be  a  complete 


Figure  2.7.  Coated  sphere  geometric  model.  Rays  are  specified  by  (p,/)  where  p  *  the 
number  of  chords  within  the  bubble,  and  /  =  the  number  of  crossings  of  the  opdc  axis. 
Dashed  lines  indicate  the  uncoated  sphere  ray  paths  while  solid  lines  represent  coated 
sphere  ray  paths.  Note  that  bending  of  the  rays  is  exaggerated  for  effect.  The  subscript  b 
indicates  rays  which  pertain  to  the  outer  sphere  only. 


description  of  the  rather  complex  phenomena  of  scattering  of  radiation  by  coated  bubbles. 
Note  also  that  the  bending  of  rays  is  exaggerated  for  illustration  of  the  effect  of  the 
coating. 

C.  Coating  Effects  on  Larger  Bubbles 

The  calculated  j  =  2  critical  angle  region  results  for  a  range  of  coating  thicknesses 
for  ka  equal  to  500, 1000,  and  2500  (a  «  37.8  pm  to  189  pm  for  =  474.6  nm) 

respectively  are  shown  in  Figs.  2.8  to  2.10.  The  parameters  used  were  similar  to  those 
used  for  Fig.  2.5.  The  effects  of  the  coating  on  increasing  bubble  size  indicate  that  as  ka 
becomes  larger  with  h  fixed  the  angular  phase  shift  in  the  coarse  structure  is  diminished. 
Mars  ton  has  calculated  this  shift  in  the  thin  film  limit  for  the  critical  angle  region  to  be 
proportional  to  h/a  which  agrees  with  the  shown  results.  Another  effect  noted  with 
increasing  bubble  size  is  that  the  changes  in  the  fine-structure  are  less  drastic  for  a  given 
value  of  h  as  bubble  size  increases.  As  an  example,  comparison  between  Figs.  2.5(e) 
and  2.9(e)  (the  3  pm  thickness  cases  for  ka  ■  100  and  ka  *  1000)  shows  that  in  the 
smaller  bubble  case  a  3  pm  thickness  (thin-dashed  line)  has  drastically  changed  the 
intensity  structure  from  the  original  noncoated  case  (solid  line)  while  for  the  large  bubble 
results,  a  lesser  shift  in  the  coarse  oscillations  and  a  slight  increase  in  fine-structure 
amplitude  has  occurred.  The  original  noncoated  structure,  however,  has  remained. 
Another  feature  overlayed  on  these  graphs  is  the  physical-optics  approximation.  The 
approximation  does  well  predicting  the  coarse  structure  for  larger  values  of  ka  for  thin 
coatings  (less  than  0.5  pm)  over  most  of  the  graphed  regions  and  is  capable  of  accounting 
for  coarse  structure  at  most  coating  thickness  values  in  the  critical  angle  region 
(approximately  80°  to  85°)  due  to  the  stationary  behavior  of  the  angular  shift  near  0C. 


Figure  2.8.  Calculated  normalized  inradiances  for  ka  -  500  (a  -  37.8  |im)  in  the  critical 
scattering  region  for  I2.  The  solid  curve  is  from  Mie  theory,  the  thin-dashed  curve  is  the 

calculated  coated  sphere  result,  and  the  thick-dashed  curve  is  from  the  phyi seal-optics 
approximation.  The  coating  thickness  values  for  *  632.8  nm  are  (a)  h  *  0.25  |im. 

(b)  h  *  0.5  nm,  (c)  1.0  pm,  (d)  h  ■  2.0  pm,  and  (e)  h  *  3.0  nm. 
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D.  Absorbing  Coatings 

Results  for  absorbing  coatings  are  shown  in  Figs.  2.1 1  -  2.14.  The  range  of 
absorbing  coatings  tested  was  n<.  =  1.5  +  ilE-9  (weakly  absorbing)  to  tig  =  1.5  +  ilE-2 

(very  strong  absorption).  The  results  in  Figs.  2.1 1  -2.13  indicate  that  the  effects  of 

absorbing  coatings  are  negligible  until  the  strong  absorbing  cases  are  reached 

(rig  =  1.5  +  ilE-3).  In  Fig.  2. 14  a  typical  result  for  the  effects  of  increasing  absorbing 

coating  thickness  is  shown.  As  can  be  seen  even  with  a  3.0  pm  coat  the  effects  are  still 
small.  This  was  even  more  apparent  for  the  perpendicular  case  which  is  given  in 
Appendix  B.  One  notable  effect  was  the  decrease  of  the  inradiance  amplitude,  which  was 
greater  at  smaller  values  of  8  and  approached  non-absorbing  results  near  the  critical 


angle.  It  should  also  be  pointed  out  that  the  coarse  structure  angular  shift  is  present  in  the 
absorbing  coating  cases  also.  Severe  angle  shift  changes  in  the  coarse  structure  amplitude 
and  structure  were  observed  in  all  cases  for  nc»1.5  +  ilE-2  (Fig.  2.13). 

In  analyzing  die  absorbing  coating  results  it  was  found  that  the  effects  of  absorption 
were  minimal.  Since  most  liquids  have  very  low  absorption  values1*  and  coating  thickness 
values  are  generally  presumed  to  be  small  it  is  doubtful  whether  absorption  will  have  a 
major  effect  on  actual  coated-bubble  measurements.  Thus,  modeling  could  be  done  using 
the  non-absorbing  cases.  This  would  speed  up  computation  time  considerably  as  the 
imaginary  components  would  not  be  involved  in  the  calculations.  Another  point  of  interest 
is  that  coatings  would  have  to  be  extremely  thick  (>  3  pm),  especially  for  larger  bubbles,  in 


order  to  have  any  major  effects  even  for  strongly  absorbing  cases.  Finally,  the  absorption 
value  for  water  is  in  the  Im(nw) *■  1 E-9  region  and  the  effects  of  such  a  small  value  of 
Im(nw)  on  the  scattering  pattern  are  expected  to  be  small.  Hence  nw  is  always  taken  to  be 
real. 


X 


Figure  2.1 1.  Calculated  normalized  irradiance  for  ka  *  500  and  =  1.5  +  ilE-3.  The 
solid  curve  is  the  non-absorbing  case  at  h  *  1.01  jim,  the  thick-dashed  curve  is  the 
physical-optics  approximation,  and  the  thin-dashed  curve  is  the  calculated  absorbing 
coating  result  for  h  -  1  01  pm. 
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(a)  h  ■  0.99  ^un,  (b)  h  ■  2.01  nm,  (c)  h  *  3.14  Jim,  and  (d)  h  a  5.00  |im.  The  solid 
curves  are  non-absorbing  results  at  the  given  values  of  h  and  the  thin-dashed  curve  is  the 
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E.  Effect  of  Varying  nc  Values 

The  fiw*!  parameter  investigated  was  the  effect  of  varying  real  values  Figures 
2.15  -  2.18  show  the  differences  between  the  =  1.5  case  generally  used  and  those  of 
tic  equal  to  1.45  and  1.55,  respectively.  In  general,  the  effect  of  coatings  with  indices  of 
refraction  less  than  1.5  was  to  diminish  the  changes  in  structure  until  the  value  n<,  =  nw  in 
which  case  the  shell  was  water  and  gave  Mie  results.  For  <  nw  an  angular  shift  of  the 
coarse  structure  in  the  opposite  direction  towards  lower  values  of  8  occurred.  In  the 
cases  where  >  1.5  the  effect  was  one  of  accelerating  both  the  shift  and  breakdown  in 
structure.  The  difference  between  values  ■  1.45, 1.5,  and  1.55  was  very  slight  for 
thin  coatings  at  all  size  parameters;  i.e.,  Figs.  2.15  and  2. 16  and  even  less  noticeable  at 
larger  values  of  ka  for  thicker  coatings  (Figs.  2.17  and  2.18). 

In  general  most  surfactants  found  in  sea  water  will  most  likely  have  indices  of 
refraction  in  the  1.45  to  1.55  range.16  Based  on  die  results  in  this  section,  the  use  of  the 
value  nw  »  1.5  will  yield  representative  results.  Although  there  are  some  small  changes 

in  the  amplitude  of  the  fine  structure,  the  amount  of  coarse  structure  angular  shift  in  this 
region  is  negligible  for  thinly  coated  bubbles  at  all  sizes  and  thickly  coated  bubbles  at  large 
values  of  ka.  Using  the  results  characterized  here  it  would  be  difficult  to  determine  from 
near-critical  angle  scattering  data  alone  if  the  coating  had  a  1.45, 1.55,  or  1.5  index  of 
refraction. 
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CHAPTER  3 


DISCUSSION 

3.1.  Interpretation  and  Analysis 

Computational  results  indicate  that  for  both  perpendicular  and  parallel 
polarizations,  I|  and  I2  respectively,  the  irradiance  patterns  of  die  scattered  light  for 

coated  bubbles  in  water  are  noticeably  altered  from  noncoated  results.  While  both  the 
coarse  and  fine  structure  oscillations  previously  evident  in  noncoated  Mie  results  remain,  a 
prominent  angular  shift  toward  8C  occurs  in  the  coarse  structure  oscillations  due  to  the 

coating.  The  degree  of  this  shift  is  dependent  upon  three  factors;  the  size  of  the  bubble  a, 
the  index  of  refraction  of  the  coating  and  finally  the  coating  thickness  h.  For  large 

bubbles  (i.e..  Fig.  2. 10)  the  amount  of  shift  is  much  less  than  for  die  smaller  sized  bubble 
(i.e..  Fig.  2.3).  The  shift  has  been  calculated  by  Professor  Marston  in  the  thin  coating 
limit  to  be  directly  proportional  to  h/a. 

Other  effects  induced  by  die  coating  are  changes  in  the  period  and  amplitude  of  the 
fine  structure  oscillations.  As  an  illustration  of  this  effect  consider  the  fine  structure 
oscillations  evident  in  Fig.  2.6  at  fixed  scattering  angle,  say  60  degrees.  In  Fig.  2.6(b) 

(h  »  0.08  4m)  the  angle  between  the  two  fine  structure  maxima  peaks  is  roughly  2.6° 
while  in  lug.  2.6(i)  (h  ■  1.12  4m)  the  distance  is  roughly  1.7°.  Thus,  in  general,  this 
indicates  that  the  fine  structure  quasi-period  A6fS  tends  to  decrease  with  increasing  h. 

The  plausible  reason  for  this  decrease  is  that  the  transverse  distance  between  the  far  side, 
or  (2,1)  ray,  and  the  surface  reflected  ray  (0,0)  appears  to  increase  with  h.  Inspection  of 
Fig.  2.6  also  indicates  that  except  for  2.6(g)  the  amplitude  of  the  fine  structure  oscillations 
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tends  to  increase  with  increasing  h.  It  is  possible  that  this  could  be  accounted  for  by  a 
shift  in  die  effective  reflection  or  transmission  coefficients  for  the  (2,1)  ray  shown  in 
Fig.  2.7.  The  observations  nwte  here  were  shown  to  be  valid  for  both  absorbing  and 
non-absorbing  coatings. 

The  computation  results  for  the  case  of  coatings  with  absorption  may  be  under¬ 
stood  as  follows:  analyzing  Figs.  2.1 1  to  2.14  one  finds  that  it  is  possible  on  physical 
grounds  to  anticipate  that  the  effects  of  absorption  are  governed  by  the  magnitude  of  a 
unitless  parameter  O  m  4fthIm(nc)Aair  The  irradiance  of  a  wave  propagating  through  an 
absorbing  medium  a  distance  h  is  related  to  0  by  the  relation  I  *  IqC"®  where  Iq  is  the 
initial  irradiance.  This  is  explained  in  Section  2.8  of  Ref.  1.  Let  us  now  illustrate  the  use 
of  this  parameter  o  by  looking  at  an  actual  case  presented  in  the  data.  Analyzing 
Fig.  2.12  we  find  that  h  ■  1.0  pm  and  ImCn^  *  5E-3  giving  a  value  of 
o  =  9.93E-2.  Calculating  a  for  Fig.  2.14(d)  with  h  *  5.0  pm  we  obtain  the  same  value 
of  o  *  9.93E-2.  Note  that  the  effect  on  the  irradiance  curves  is  similar.  Finally,  it  should 
be  noted  that  for  higher  o  values  the  effects  due  to  absorption  become  more  drastic. 
Comparing  Fig.  2.13  (o  *  1.99E-1)  with  Figs.  2.11  and  2.14(a)  (a  *  1.99E-2)  we  find 
this  is  indeed  the  case.  Thus  for  coatings  with  0  less  than  approximately  0.02  the  effects 
of  absorption  are  negligible  and  real  values  for  would  best  be  used  in  computations 

since  the  use  of  real  values  reduces  the  computation  time  significantly. 

Based  on  the  results  presented  here  and  in  Chapter  2,  the  following  conclusions 

/ 

are  noted  for  the  case  of  j  *  2  scattering  from  bubbles  having  nonabsorbing  coatings. 
First,  for  a  coating  less  titan  1  pm  in  thickness  on  bubbles  having  ka  2  500  there  is  very 
little  effect  of  the  coating  on  the  principal  features  of  the  scattering  pattern.  As  noted  in 
Appendix  B,  the  effect  of  a  given  coating  on  the  coarse  features  of  the  scattering  pattern  is 
usually  smaller  for  j  ■  2  than  for  j  =  1 .  This  may  be  because  the  reflection  of  light  from 
the  coating-water  interface  is  smaller  for  j  *  2  than  j  =  1  in  the  angular  region  of 


interest  (See  Appendix  B.)  For  bubbles  with  ka  ■  100,  changes  in  the  fine  structure  can 
be  pronounced  for  thickness  values  between  0.25  |im  and  1.0  nm. 


3.2.  Experimental  Applications 

If  microbubbles  in  the  ocean  are  proven  to  be  coated  by  a  thin  film  of  surfactant 
then  it  becomes  important  to  determine  if  such  parameters  as  thickness  and  complex 
refractive  index  erf  the  coatings  affect  the  reliability  where  by  light  scattering  techniques 
can  be  used  to  determine  die  radius  a  of  a  gas  bubble.  Calculations  given  here  for  both 
coated  and  noncoated  bubbles  indicate  that  for  their  sizing  and  detection  the  coarse 
structure  oscillations  for  the  j  =  2  case  in  the  critical  angle  scattering  region  are  best  used 
for  sizing.  The  basis  for  this  observation  is  three-fold;  first,  the  fine  structure  is  generally 
much  smaller  in  magnitude  for  both  coated  and  noncoated  bubbles  for  j  *  2  than  for 
j  =  1.  Secondly,  the  j  *  2  case  for  coated  values  does  not  have  the  change  in  coarse 
structure  (other  than  the  angular  shift  which  was  noted  for  both  cases)  as  was  observed 
for  j »  1;  see  e.g..  Fig.  B.3(e).  Thirdly,  unless  coatings  are  found  to  be  quite  thick  or 
highly  absorbing  previous  observations  and  conclusions  made  using  Mie  theory  and  the 
physical-optics  approximation  could  be  used  as  an  aid  in  the  general  design  of  apparatus. 
Thus,  it  is  recommended  that  sizing  and  detecting  of  coated  bubbles  in  experimental 
applications  is  best  done  by  making  use  of  the  j  *  2  coarse-structure  and  placing  the 
maximum  number  of  detectors  as  close  as  possible  to  the  critical  angle  scattering  region. 


APPENDIX  A 


COMPUTER  PROGRAMS 

Two  programs  were  used  to  generate  our  plots.  COATSPHR,  using  modified 
versions  of  Wiscombe's  Mievo,  and  Bohren  and  Huffman's  coated  sphere  subroutines 
produced  the  coated  sphere  results  and  ESCAT  the  approximation  results.  The  data  was 
stored  in  system  files  and  manipulated  using  a  Statistical  Analysis  System  (SAS)  routine 
which  generated  the  plots  presented.  All  plots  were  done  on  an  HP7475a  plotter  using 
0.5  mm  and  0.7  mm  Hewlett  Packard  drafting  plotter  pens.  All  programs  are  written  in 
IBM  FORTRANVS. 

COATSPHR 

To  use  COATSPHR  the  following  values  must  be  input  in  the  main  calling 
program: 

1 .  DT  is  the  angle  increment  to  be  sampled.  Typically  it  must  be  less  than  (A/7a) 
rad  where  a  is  the  inner  sphere  radius. 

2.  ANGfN)  is  the  initial  angle  at  which  scattering  calculations  begin.  All 
calculations  are  symmetric  about  90“  meaning  that  if  ANG(N)  *  30°  the  range 
of  30°  to  150°  will  be  calculated. 

3.  NANG  is  the  number  of  angles  for  which  the  coefficients  are  calculated. 
Typically  for  the  entire  range  it  is  determined  by  NANG  =  1/DT  x  (90°  - 
ANG(N))  +  1.  Note  that  for  this  particular  program  version  NANG  must  be 
less  than  1800.  (This  number  may  be  increased  by  increasing  all  array  sizes 
but  for  the  cases  presented  1800  was  adequate.)  If  only  a  small  region  is 


desired  then  the  variable  NANG  should  be  determined  by  NANG  =  1/DT  x 
Cast  angle)  +  1  (for  example,  to  compute  0°  to  5°  and  175°  to  180°  with 
DT  -  1/20  use  NANG  -  20  x  5  +  1  -  101;  ANG(N)  -  0°. 

In  addition  to  this  the  following  values  must  be  input  in  subroutine  BHCP(X,Y). 

1 .  REFMED  is  the  refractive  index  of  the  surrounding  medium. 

2.  REFRE1  is  the  real  part  of  the  index  of  refraction  of  the  inner  sphere  medium. 

3.  REFTM1  is  the  imaginary  component  of  the  index  of  refraction  of  the  inner 
sphere  medium. 

4.  REFRE2  is  the  real  part  of  index  of  refraction  of  the  coating. 

3.  REFTM2  is  the  imaginary  component  of  the  coating  index  of  refraction. 

6.  AOB  is  the  ratio  of  the  inner  radius  a  divided  by  the  outer  radius  b  (see  Fig. 
1.1).  The  coating  value  is  given  by  the  relation  h  =  b(l  -  AOB). 

7.  X  is  the  inner  sphere  size  parameter  ka  »  Ito/Kq  where  Xq  is  the 

wavelength  in  the  outer  medium. 

Once  these  values  are  input  COATSPHR  generates  the  following  output: 

1.  Y  is  size  parameter  of  the  outer  shell  of  radius  b. 

2.  Ij(0)  is  the  intensity  for  the  electric  field  perpendicular  to  the  scattering  plane. 
It  is  given  for  each  value  of  0  input  It  is  normalized  such  that  lj(0)  =  1 

represents  perfect  reflection  according  to  geometric  optics. 

3.  12(6)  is  the  intensity  for  the  parallel  case. 

4.  S|(0)  and  S2(6)  are  the  partial  wave  angular  scattering  functions. 

3.  The  average  intensity. 

6.  Efficiency  factors: 

a)  Extinction  coefficient  Qext  =*  C^/G  where  is  the  extinction  cross 
section  and  G  is  defined  below. 


b)  Scattering  coefficient  Qjc*  »  where  =  scattering  cross  section 

and  G  is  the  panicle  cross  sectional  area  projected  onto  a  plane 
perpendicular  to  the  incident  beam  (e.g.,  G  3  jtb 2  for  coated  sphere). 

c)  Absorption  coefficient  which  is  simply  0CTt  -  Qy)l. 

7.  The  asymmetry  factor  is  the  average  cosine  of  the  scattering  angle.  Fora 
particle  that  scatters  light  isotropically  the  asymmetry  factor  would  be  zero. 
Being  that  the  author  is  a  novice  in  FORTRAN  programming  no  claim  is  made  as 
to  the  efficiency  or  lack  thereof  in  COATSPHR.  The  program  requires  close  to  1 
Megabyte  of  virtual  memory  to  run  as  there  are  3  large  arrays  and  the  use  of  quadruple 
precision  to  contend  with.  Fortunately  this  is  well  within  the  IBM  3090-200's  capability 
range  and  so  no  attempt  was  made  to  economize  the  program. 

ESC  AT 

Values  for  the  physical-optics  approximation  were  generated  using  ESCAT 
FORTRAN.  The  program  is  virtually  unchanged  from  the  version  used  by  Langley  in 
1982.  Only  a  few  modifications  were  necessary  to  convert  it  for  use  on  the  FORTVS 
compiler.  To  use  ESCAT  the  following  values  must  be  entered: 

1.  M  is  the  relative  refractive  index  of  the  non-coated  bubble. 

2.  KA  is  size  parameter  of  the  noncoated  bubble. 

3.  PHEBGN  is  the  value  of  the  first  scattering  angle  (in  degrees)  to  be  evaluated. 

4.  PHIEND  is  the  last  scattering  angle  (in  degrees). 

3.  DPHI  is  the  angle  step  size  (analagous  to  DT  in  COATSPHR). 

The  output  gives  the  irradiances  Ij  for  both  the  j  *  l  and  j  ■  2  cases,  with  or 

without  diffraction  effects.  The  program  only  determines  the  coarse  structure  as  has  been 
shown  in  the  figures  but  has  the  advantage  that  it  executes  extremely  quickly  for  all  bubble 
sizes,  whereas  similar  Mie  calculations  require  significantly  more  computation  time. 


non  n  ooooo  o  n  nnnnn  nnnnnn 


FILE:  COATS PHR  FORTRAN  A. 


SCATTERING  PROGRAM  FOR  COATSO  SPHERES 
3.  C.  BILLETTS  NOVEMBER  ISIS 
ANGULAR  DIMENSIONS 
COMPLEX**  SB ACT.  SI.  S2 

REAL* IS  ANG (3600)  , XMU , PI , CON , DT , XX , CQNV ,  YY 
REAL*1S  II. 12, INTSN, ANGLE. N2 CUT 

COMMON/INOUT/S 1(3600)  ,32(3600)  ,XMU(3600)  ,  SBACE ,  XX ,  YY , 

•  QEXT , Q3CA , GQS  C , NUMANG , X2CUT 

USE  SUBROUTINE  3KC7  TO  INITIALIZE  FOR  AN  AND  BN  CALCULATIONS 
AND  OBTAIN  CORE  AND  COAT  SIZE  PARAMETERS 

CALL  SHCP(XX.YY) 

N2CUT  •  O.OQO 

PI  -  3 . 141392633 S8979323S462S433S3279SQ0 

con  •  px/iao.ao 

NANG  -  301 
CONV  -  XX**2 
DT  «  1.0QO/20.0QO 
N  -  1 

ANG(M)  *  T5.0Q0 

00  1  J  -  1,  NANG 

XMU (N)  -  QCOS(CON*ANG(N) ) 

ANG<N-*1)  »  ANG(N)*OT 

1  N  •  N+l 
NUMANG  -  2«(N-1) 

J  -  NUMANG-*  1 

DO  2  I  -  N, NUMANG 

ANG(I)  -  180 . -ANGfJ-I) 

2  XMU(t)  -  -XMU(J-I) 

CALL  SUBROUNTINE  MIS70  TO  DETERMINE  SI  AND  52 
CALL  MZEVO 
CONTINUE 

CABS  -  QEXT-43CA 
WRITE (6, 1000)  xx, a 

00  10  I  -  1, NUMANG 

11  ■  4. OQC* ( : REAL (SI (I) ) ) ••2*(AIMAG(S1 'I) ) ) **2 ) /CONV 

12  •  4.0QOV :RXAL;S2(I) ) ) ••2-»(AXMAG(S2(l; ) ) **2)/C0NV 
INTER  -  .3Q0«(I1~I2) 

ANGLE  -  ANG ( I ) 

WRITE  STATEMENTS 


10  WRITE  (6 ,  1001)  ANGLE. II. 12, SKI)  ,S2(I)  ,  INTEN 
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WRITS(4,ia02)  QEXT,QSCA,QABS,S«SC 

WRITS (4.1003)  SBACX 

STOP 


FORMAT  STATEMENTS 

.000  FORMAT ( 5X, 'CORE  PARAMETSR  — ' , D14. 4 , 3X. ' COAT  PARAMZTSR  -',314.4, 

•  //,'  ANGLZ' , SX, 'll', 11X, ' 12 ' , 20X, 1  SI' , 24X, 

•  '32 ' , 1SX, ' INTENSITY' ) 

.001  FORMAT (FS . 3 , 8E14 .4) 

.002  FORMAT (  /30X,  'EXTINCTION  SCATTERING  ABSORPTION ’/?X, 

•  ' EFFICIENCY  FACTORS ', 3F14 . 4/  ?X, 'ASYMMETRY  FACTOR  -\F9.4) 
.003  F0RMAT(7X,'S  BACKS CATTSS  -'.2E14.4) 

END 

SUBROUTINE  MXEV0 


NOTE:  THIS  SUBROUTINE  USES  THE  SAME  VARIABLE  NAMES 
AS  THOSE  FOUND  IN  WISCOMBE'S  MXEV0  FROM  WHICH  IT 
IS  DERIVED.  SEX  MIEVO  FOR  A  COMPLETE  LISTING. 


COMMON/INOUT/Sl ( 3400) ,S2 (3400) , XMU(3400) , SBACK, XX, YY, 

•  QEXT , QSCA , SG3C , NUMANG , N2CUT/ BLOCE2/AN , BN/ BLOCX3/NT 
RIAL* 14  NP1DN, N2CUT 

REAL*14  XX, YY, XINV , FN, FNP1 , RN, RNP1, TWONP1 , COEFF 

REAL* 14  PIN. PZNM1. IMP. TACK, XKU 

LOGICAL  NOIMAG,  HOANGS 

COMPLEX* 4  3BACX.  SI,  S3 

COMPLEX* 3 2  SP.SM.SPS.SMS 

COMPLEX* 3 2  ANPM.BNPM,ANP,BNP,ANM1.SNM1 

COMPLEX* 3 2  AN, 3N 

DIMENSION  SP(ISOO) ,SM(1400) ,S?S(1800) ,SMS(1I00) , 

*  PIN ( 1400) , RINM1 ( 1400) ,TAUN(1800) ,TMP(1S00) 

XINV  -  l.OQO/XX 
NN2  -  NUMANG* 1 
NN  -  NN2/2 

HOANGS  -  nCtMANG  .  ZQ .  0 

CALCULATE  NUMBER  OF  TERMS  IN  PARTIAL  WAVE  SERIES 
USING  EMPIRICAL  FORMULAE  WHICH  WERE  FITTED  FOR 
SIZE  PARAMETERS  UP  TO  20,000 

IF ( YY . LE . I . 0Q0 )  NT  -  YY*4 . 00Q0*YY** ( 1 . OQO/3 . OQO) -1 . 0Q0 
1F( YY.GT. 4 . 0Q0)  NT  •  YY*4. 03Q0*YY** ( 1. 3Q0/3 . OQO) -2 . 0Q0 
HTP1  -  NT*i 
HTP2  -  NT *2 

WRITE (4,4001)  NT 
1001  FORMAT (4X,  '  NT  -M4/4X) 


INITIALIZE  QUANTITIES  USEO  TOR  EFFICIENT  CALCULATION  OF 
NUMERICAL  COEFFICIENTS  IN  PARTIAL  WAVE  SERIES 
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AN PM  -  MM*ANP 
3NPM  -  MM* BNP 

C 

00  400  J  -  1,NN 

ADO  OP  SUMS  WHILE  UPWARD  RECURS INS  ANGULAR  rUNCTIONS 
LITTLE  ?I  AND  LITTLE  TAU 

TMP(J)  -  (XMU C J ) *PIN( J) )  -  PINMl(J) 

TAUN(J)  -  FN*TMP(J)  -  PINMl(J) 

SP(J)  -  SP(J)  *  ANP*(PIN(J)-TAUN(J) ) 

SMS (J)  -  SMS (I)  *  3NPM*(PIN(J)-TAUN(J) ) 

SM(J)  •  SM(J)  *  3NP*(PIN(J)-TACH(J) ) 

SPS(J)  -  SPS(J)  *  ANPM*(PIN(J) -TAUN(J) ) 

PINMl(J)  -  PIN{ J) 

PIN (J)  -  (XMU(J) *PIN(J) )  *  NP1DN*TMP(J) 

400  CONTINUE 

UPDATE  RELEVANT  QUANTITIES  FOR  NEXT  PASS  THROUGH  LOOP 

4S0  MM  -  -MM 
FN  •  FNP1 
RN  -  RNP1 
ANM1  -  AN 
3MK1  -  3N 
500  CONTINUE 

MULTIPLY  SUMS  3Y  APPROPRIATE  FACTORS  TO  GET  QEXT,  QSCA,  GQSC 

QEXT  -  2 . 0Q0* (XINV**2 ) *QEXT 
QSCA  -  2 . 0Q0* ( XINV**2 ) *QSCA 
IF(NOIMAG)  QSCA  -  QEXT 
3QSC  •  4 . 0Q0  * ( XINV *  «  2 ) *GQSC • 

SHACK  »  0.3Q0*SBACK 

IF (N0ANG3)  RETURN 
30  800  J  »  I, NM 

S1(J)  w  0 . SQO* (3P( J) »SM( J) ) 

S2(J)  -  0 . SQO* (SP( J) -3M( J) ) 

SICNN2-J)-  a.SQO^SPSdl-SMSCJ)  ) 

BOO  32 (NN2-J) •  0.3Q0*(S?S(J) -SMS(J) ) 

RETURN 


END 

SUBROUTINE  BHCP ( X , Y ) 

SUBROUTINE  TO  INPUT  COATED  SPHERE  PARAMETERS  AND  TO 
INITIALIZE  FOR  LATER  ENTRY  INTO  BC0AT1 

REAL* 18  X,  Y, PI , AOB 
REAL* IS  QEXT, QSCA, Q8ACX 

REAL* IS  REFMEO , REFREI , REFXM1 , REFRE2 , REFIM2 
COMPLEX*] 2  RFREL1,  RFREL2 

REFMEO  -  REFRACTIVE  INDEX  OF  SURROUNDING  MEDIUM 
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refractive  otex  or  core  -  refrei  -  :*refi>ii 

REFRACTIVE  INDEX  OF  COAT  -  REFREl  -  I* RE FIN 2 
A08  -  RATIO  OF  CORE  RADIUS  TO  COAT  RADIUS 
X  -  CORE  SIZE  PARAMETER 
¥  -  COAT  SIZE  PARAMETER 

REFMED  -  4 . OQO/3 . 0Q0 
RETREI  -  1.00Q0 
RIFIM1  -  O.OOQO 
REFRE2  -  I.SOQO 
REFIM2  -  O.OOQO 
AOB  -  0.99Q0 
X  -  100. OQO 
¥  -  X/AOB 

DEFINE  VARIABLES  RFREL1 ,  RTREL2 

RFREL1  -  QCM7LX(R£FRE1,REFIM1) /REFMED 
RFREL2  -  QCMPLX ( REFRE2 , REFIM2 ) / REFMED 

USE  SUBROUTINE  BHCOAT  TO  CALC.  EFFICIENCIES  FOR  COATED  SPHERES 

CALL  BHCOAT  (X, Y, RFR1L1 , RFREL2 ) 

WRITE  STATEMENTS 

WRITE  (4,9) 

WRITE  (4,11) 

WRITE  (4,9) 

WRITE  (4,4) 

WRITE  (4,12)  REFMED , REFREl , RZFTM1 , RZFRE2 , REFIM2 
WRITE  (4,14)  AOB 

FORMAT  STATEMENTS 


4  FORMAT ('  ') 

9  F0RMAT'19X, **••••*••*••••••••••••*••••***•**• ) 

11  FORMAT (19X, 'COATED  SPHERE  SCATTERING  PROGRAM* ) 

12  FORMAT  (SX. ’REFMED  -*,  E14 . 9/ SX,  *  REFREl  -\S14.9, 

•JX.'REFIMl  •' , S14 . 9/SX, 'REFRS2  ■' , £14. 9, OX, 'REFIM2  -*,214.9) 

14  FORMAT (SX, 'AOB  -*,E14.9) 

RETURN 
97  3TOP 

END 

SUBROUTINE  BHCOAT  ( X , ¥ , RFREL1 , RFREL2 ) 

SUBROUTINE  BHCOAT  TO  CALCULATE  LITTLE  AN  AND  BN  FOR 
GIVEN  SIZE  PARAMETERS  OF  CORE  AND  COAT  ANO  RELATIVE 
REFRACTIVE  INDICES.  BASED  ON  SUBROUTINE  GIVEN  IN  AP¬ 
PENDIX  B  OF  ABSORPTION  ANO  SCATTERING  OF  LIGHT  B¥  C. 
BOHREN  ANO  0.  HUFFMAN  (WILEY,  19(3) 

*••••••••*•••••• •••••••CAUTION****** ••••*•*•*••*••••• 

BHCOAT  SHOULD  NOT  BE  USED  FOR  LARGE,  HIGHLY  ABSORBING 


E:  CCATSPHR  FORTRAN  A1 


COATS D  SPHERES  X*R£EIM1,  X»R£EIM2  AND  Y*R£7IM2  should 
32  L2SS  THAN  ABOUT  30. 

. . 

COMMON/BLOCX2/AN, 3N/BLOCK3/YSTOP 
REAL * 1 8  QSCA.Q BACK , 3 EXT, PSIY , CHIY , DEL, SN, YSTOP 
REAL*1S  X,Y,?SI0Y,?SI1Y,CH23Y,CHI1Y 
COMPLEX*! 2  RESELL , REREL2 , XL , X2 ,  Y2 , REFREL 
COMPLEX*32  01X1 , 30X1,01X2 , 30X2 , D1Y2 , 30Y2 

COMPLEX*!!  XT3Y ,  XII Y  , XZY, CHI0Y2 , CHI1Y2 , CHIY2 , CHX0X2 , CHI1X2 , CHTX2 
COMPLEX* 3 2  CHIPX2 , CHI? Y2 , AN CAP , BNCAP , 3NBAR , GNBAR , AN , BN , CRACK , 3 RACK 
COMPLEX* 3 2  XBACK.AMESS1, AMESS2 , AMESS3 ,AKESS4 . 

COMPLEX* 3 2  CBICOl , 8IG01, 0BI5D2 , 3IGD2 , 3SIGD3 , 3IGD3 
COMPLEX* 3 2  ZINV, ZINV2 , 2INV3 

DIMENSION  BIGOI(IOISO) ,3IGD2(10130) , BIG03 (10150) 

DEL  IS  THE  INNER  SPHERE  DIVERGENCE  CRITERION 

DEL  -  1.0q-8 
XI  *  RERELL*X 
X2  -  RFREL2*X 
Y2  -  RTREL2*Y 

YSTOP  -  Y  *  4 . 35Q0*Y** ( 1 . OQO/ 3 . OQO)  -  2.0QO 

REFREL  *=  REREL2  /  RFREL1 

NSTOP  -  YSTOP 

NTP1  •  NSTOP  *  1 

NTP2  -  NSTOP  -  2 

SERIES  TERMINATED  AFTER  NSTOP  TERMS 

SUBROUTINE  TO  CALCULATE  INITIAL  HIGH  OROER  3IG01-N 

CALL  3IG0 (X, RESELL, OBIGOl) 

3IG01 (NSTOP)  -  OBIGOl 
ZINV  -  1. OQO/ (RESELL *X) 

30  70  N  ■  2, NSTOP 

OBIGOl  -  (NTP2-N) *ZINV-1. OQO/ (( (NTP2-N) *ZINV) -OBIGOl) 

*3  3IG0KNTP1-N)  -  OBIGOl 

SUBROUTINE  TO  CALCULATE  INITIAL  HIGH  OROER  3IGD2-N 

CALL  BXGO(X, RER2L2 , 0BIG02 ) 

BIG02 (NSTOP)  •  0BIG02 
ZINV2  -  l.OQO/ (SFR£L2*X) 

00  «0  N  -  a, NSTOP 

3BIGD2  -  (NTP2-N) *ZINV2-1 . OQO/ ( ( (NTP2-N) *ZINV2.) -DBIGD2) 

80  BIG02 (NTPl-N)  -  081002 


SUBROUTINE  TO  CALCULATE  INITIAL  HIGH  OROER  3IG03-N 
CALL  SIGO(Y, RESELL . 0BIG33 ) 


96 
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SI003(HSTOP)  -  38IC03 
2INV3  -  l.OQO/  (RFREL2*Y) 

DO  90  K  ■  3 , NSTOP 

08T003  -  (NTP2-N) *IXNV3-1.0Q0/( ( (NTP2-N) *TZNV3)-08t003) 
90  BI0D3 (NTP1-N)  -  38X003 

00X1  -  CQCOS (XI) /CQSIN (XI)  * 

00X3  ■  CQCOS ( X2 )/CQS IN (X2) 

00Y2  -  CQCOS ( Y2 )/CQ3 IN (Y2) 

PSIOY  -  QCOS(Y) 

PSX1Y  -  QSIN(Y) 

CHXOY  -  -QSXN(Y) 

CHX1Y  -  QCOS(Y) 

XXOY  -  QCMFLX( PSIOY, -CHXOY) 

XI 1Y  -  QCMPLX(PSI1Y,  -CHI1Y) 

CHZ0Y2  -  -CQSIN  ('.'2 ) 

CHI1Y2  -  CQCOS (Y2) 

CHI 0X2  -  -CQSIN (X2) 

CHI1X3  -  CQCOS (X2) 

QSCA  •  O.OQO 
QEX7  -  O.OQC 

XBACX  -  Q CMP LX  (0 . OQO , 0 . OQO) 

N  «  1 
I FLAG  -  0 

RETURN 
ENTRY  BC0AT1 
200  SN  -  N 

PSIY  -  (2.0Q0*RN-1.0Q0) ‘P3Z1Y/Y  -  PSIOY 
CHIY  »  ( 2 . OQO *RN-1. OQO) *CHI1Y/Y  -  CHXOY 
XIY  «  QCMPLX( PSIY, -CHIY) 

01Y2  -  81003(H) 

IF  (IFLAC.EQ.l)  00  TO  999 
31X1  ■  BtCDl(N) 

01X2  •  BZ0S2(N) 

CHIX3  •  ( 2 . OQO*RN  -  l.OQO)*  CHI1X2/X3  -CHI 3X2 
CHIY 2  -  (2 . OQO*RN  -  l.OQO)*  CHI1Y2/Y2  -CHI0Y2 
CHI 7X2  •  CHI 1X2  -  RN*CHIX2/X2 
CHIPY2  •  CHI1Y2  -  RN*CHIY2/Y2 

AN CAP  -  RXTRIL* 01X1-01X2 

ANCAP  ■  ANCAP/ (RIFREL*D1X1*CHIX2  -  CHXPX2 ) 

AN CAP  -  ANCAP/ (CHIX2 *01X2  -  CHIPX2) 

BRACK  -  ANCAP* (CHIY2*01Y2  -  CHIPY2) 

BNCAP  -  RZFREL*01X2  -  01X1 

BN CAP  ■  BNCAP/ (RXFRXL*CHI?X2  -  01X1*CHIX2) 

BNCAP  •  BNCAP/ (CHIX2 *01X2  -  CHIPX2) 

CRACK  •  BNCAP* (CHIY2*01Y2  -  CHIPY2) 


it  Jt 


■l  JO. 


FILE:  COATS? SR  FORTRAN  A* 


AMESS1  -  3RACX*CHI?Y2 
AMESS2  -  3RACX*CHIY2 
AKESS3  -  CRACX*CHI?Y2 
AMISS 4  -  CRACK* CHI Y 2 

IF  (C0A8S  (AMESS1)-.  ST.  QEL*C0ABS  ( D1Y2) )  SO  70  999 
I? (COABS ( AME5S2 ) . ST.2EL)  SO  TO  999 
I? (COABS (AKE5S3 ) . ST. 2EL*CQABS (01Y2) )  SO  TO  999 
IF(CQABS(AKESS4) .ST.2EL)  SO  TO  999 

BRACK  -  QCMPLX  (0.300,9.900) 

CRACK  -  QCMPLX  (0.900,3.900) 

I  TEAS  -  1 

999  0N1AK  -  0112  -  BRACK*CHt?Y2 

ONBAR  -  ONBAR/ (1. 9Q0  -  3FACX*CHIY2) 

SNBAR  -  01Y2  -  CRACK»CHXPY2 
SHEAR  -  SHEAR/ (1.000  -  CRACK*CHIY2) 

AM  -  (DHBAR/RFREL2fRH/Y) *PSIY  -  PSI1Y 
AH  •  AM/( (OHBAR/RFRIL2+RH/Y) »XIY  -  XX1Y) 

3M  -  (RFREL2*SHBAR»RH/Y) *PSIY  -  pSIIY 
3H  -  3M/( (RTR£L2«GNBAR+RN/Y) *XIY-XI1Y) 

PSIOY  -  PSX1Y 
PSI1Y  -  PSXY 
CHIOY  •  CHI1Y 
CHI1Y  -  CHIY 

XX 1Y  -  QCMPLX ( PSIIY. -CHX1Y) 

CHI 9X2  -  CHI 1X2 
CHI 1X3  -  CHIX2 
CHI3Y2  -  CHI1Y2 
CHX1Y2  »  CHIY2 

20X1  ■  21X1 
20X2  -  21X2 
D0Y2  -  01Y2 
M  •  H  -  1 

RETURN 

SHO 

SUBROUTINE  BIS0(XX, XOR, OBISO) 

SUBROUTINE  TO  CALCULATE  LOGARITHMIC  2ERIVATIVSS  3Y  LZNZ’S 
KETHOO  US INS  20WN  RECURSION 

CCMMOH/ BLOCK3 / YSTO  P 
REAL* 14  XINV.XX, YSTOP 

COMPLEX* 3 2  ZINV , IOR , FF , AX , SEN , HUM , NTH , DTD , TT , 28XSD 
OATA  EPS1/1.Q-2/EPS2/ 1.0-9/ 

DATA  MAXIT/19000/ 

NT  •  YSTOP 
HTP1  -  NT  -  1 
NTP2  •  NT  *  2 


ono  n  o  onoo  ooo  oo 
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FILE:  COATS PHR  FORTRAN  A1 


MAXE  StJR£  9XS0  ARRAY  WILL  BE  SIS  ENOUGH 

ITCNTPl.LE. 20130)  GO  TO  10 
WRITE («,«000)  NT 

1000  FORMAT (///'ESTIMATED  LENGTH  Or  HIE  SERIES  NT-', IS) 

STOP 

COMPOTE  SISO 

10  XI NV  >  1.0Q0/ XX 
ZINV  -  XINV/IOR 

PREPARE  FOR  DOWN  RE  CURRANCE •COMPUTE  INITIAL  HIGH  ORDER 
8IG0-N  USING  LENTZ  METHOD 

FT  -  NTP1*ZINV 

MM  -  -1 

30C  -  2*NT*3 

AX  -  (MM* XX) *ZINV 

DEN  -  AX 

NCM  -  OEN  *  1.0QO/FT 
FOUNT  -  1 

20  KOUNT  -  FOUNT*! 

IT f FOUNT . GT . MAXIT)  GO  TO  40 

IT(CQAB3 (NUM/AX) . GT . EPSI . AND. CQABS ( DEN/AK) .ST.SPSI)  GO  TO  30 

MM  -  -MM 

XX  •  KK*2 

AX  -  (MM*KX)*ZINV 

MTN  -  AK*NCM  -  l.OQO 

DTD  -  AX* OEN  -  l.OQO 

FT  -  (NTN/DTD) *FF 

MM  -  -MM 

XX  -  KX*2 

AX  -  (MH*KK) *ZINV 

HUM  -  AX  *  NUM/NTN 

DEN  -  AX  -  OEN/ DTD 

FOUNT  -  KOUNT  *  1 

GO  TO  20 

30  TT  -  MUM/ OEN 
FT  •  TT*FF 

CHECX  FOR  CONVERGENCE 

ET(QASS(QRXAL(TT)-t.OQO) .LT.SPS2  .AND.  QASS (QIMAC (TT) ) . LT. £?S2 
•  GO  TO  SO 
MM  •  -MM 
XX  -  XX*  2 
AX  -  (MM»KX) *ZIMV 
NUM  ■  AX  -  1. OQO/NUM 
OEN  -  AX  -  1.0Q0/DEN 
SO  TO  20 

40  WRITE(S , (001)  NT, XX, I0R, AX/ NUM. DEN, TT, FT 


TIX:  COATSPHR  FORTRAN  A1 


3001  FORMA? (///'CONTINUED  FRACTION  FOR  A-SUB-NT  FAILED  ?3  CONVERSE ' 

•  '  NT«',I6/’  X-* ,£20.3/ '  RIFR  INDEX- ' , 2E20 . 3/ '  AK-',2S20.3/ 

•  '  MUM- 1 , 2E20 . 3/ '  DEM-1 , 2S20 . 3/ *  IT-' , 2E20.3/ ' FF-‘ , 2E23.3) 
STOP 

so  obisd  -  n 

RETURN 


OUUUOtJUUUUUUUUOUOOU  UUU  «>  If  u  uuu  uou  o  u 


FILE:  ISCV 


FORTRAN  A. 


PROGRAM  S3 CAT  TO  COMPUTE  COMPLEX  SCATTERING 
AMPLITUDES  SI  AMO  S3  USIMG  A  PHYSICAL  OPTICS 
APPROXIMATION  AND  DIFFRACTION. 

O.L.  KINGSBURY  AND  P.L.  MASSTON  1910 
O.S.  LANGLEY  19S3 

MOOIITED  FOR  CMS  BY  S.C.SILLZTTE  APRIL  7,  1986 


INPUT  PARAMETERS 


M  -  RELATIVE  REFRACTIVE  INDEX  OF  BUBBLE 
KA  -  SHE  PARAMETER 
PRISON  -  STARTING  ANGLE  IN  DECREES 
PKISNO  -  LAST  ANGLE  (DEGREES) 

DPMI  -  ANGLE  STS?  SIZE  (DEGREES) 


REAL  M, KA, II , 12 , ID , IT1  IT2 

COMPLEX  si ,  Si', sal,  S02 Isll ,  S12 ,  F01 ,  F02 .  ru ,  F12  .  SOI ,  SD2 
DOUBLE  PRECISION  N.TRC.ST1 
DATA  RMIN/0.0/ 

H  -  .78 
KA  -  100.0 
PRISON  >30.0 
PHI END  >90.0 
DPHI  >  1.0/18.0 


II  >  1  GIVES  P.O.A  RESULTS  USING  ONLY  (0,3)  i  (1,0)  RAYS 

II  •  2 

UJ  >  I  GIVES  DIFFRACTION  ONLY  RESULTS 

TU  •  0 

XX  >  1  GIVES  P.O.A.  INCLUDING  DIFFRACTION  RESULTS 

KX  >  0 

WRITS (1.1000)  K, KA 

1300  FORMAT ('M  •  1 ,F6.4,'  KA>’,F6.2) 

COMPUTE  CRITICAL  ANGLES  FOR  GIVEN  REFRACTIVE  INDEX 

PI  •  4 . 000*DATAN f 1 . 300 ) 

N  -  1.000 /M 

THC  ■  DATAN ( 1 . 3D0/DSQRT (N**2  -  1.000)) 

PRC  >  PI  -  2 . 0*THC 

BEGIN  COMPUTATION  FOR  SCATTERING  ANGLE  PKI  (RADIANS) 


L  •  0 

PRIO  ■  PRISON 
IF ( PRID . ST . PKISNO :  STOP 
30  PHI  >  PHID*PI/i30.: 


COMPUTE  ( 3 . :  RAY  PARAMETERS 


/•  .%  .  / 


nnoonnnnn  n  o  n  n  o  o  n  n  n  o  n  n 


FTLZ:  SSCAT 


FORTRAN  A1 


tho  -  incidence  angle 

501,502  -  PHASES  FOR  J  -  1.2  ; CENTER  RET.) 

F01.F02  -  S.P.A.  FAR  ZONE  COMPLEX  AMPLITUDES  FOR  J  -  1,2 

THO  -  (PI-PHI) /2. 3 
CTO  -  COS (THO) 

STO  •  SIN (THO) 

W  -  SIN (PHC-PHI) *SQRT  ( JCA*COS ( SNSL(THC) ) / ( 2 . 0*PX) ) 

CALL  FREE (W, FC, FS) 

F01  -  CMPLX(FC.FS) *CSX? (CMPLXfO . 0 , -1 . 3*PI/4 . 0) )/SQRT(2.0) 

F02  -  FOl 

GOl  -  -2 . 0»KA»CT0 

502  -  501 

INITIALIZE  (1.0)  RAY  SCATTERING  AMPLITUDES: 

IP  PHI  EXCEEDS  PHC  THERE  IS  NO  (1,0)  RAY,  AND  NO 
(0,0)  REFLECTION  PRASE  SHIFT. 

511  -  CMPLX(0. 0,0.0) 

512  -  CMPLX(0. 0,0.0) 

IF(PHI.5E.PHC)  SO  TO  SO 

COMPUTE  (1.0)  RAY  PARAMETERS: 

TH1  -  INCIDENT  ANGLE 
ROl  -  REFRACTION  ANGLE 

GU.512  -  PHASES  FOR  J  -  1,2  (CENTER  REFERENCED) 
D1  -  GEOMETRIC  DIVERGENCE  FACTOR 
R1I.R12  -  FRESNEL  COEFFICIENTS  FOR  J  •  1.2 
FU.F12  -  S.P.A.  FAR  ZONE  COMPLEX  AMPLITUDES 

PH2  -  PHI/2.0 

TH1  •  ATAN (M*3IN(PH2 ) / ( I. 3-M*C0S ( PH2) ) ) 

CT1  -  COS(THl) 

ST1  -  DSIN (DBLE (TH1) ) 

R01  -  OATAN(N*ST1/DSQRT(1.000-(N*STL) **2) ) 

CR1  -  COS (R01) 

511  -  2.0*KA*(M*CR1-CT1) 

512  -  511 
TAU  -  N 

IF ( CR1 . 5T . RMIN)  TAU  -  TAU*  CT1/CS1 

01  •  K*CT1**2/ ( * . 0* (TAU-1. 0 ) • ( 1 . 3-M*C0S ( PH2 ) ) *C0S ( ?H2 ) ) 

0IF1  -  ROl  -  TH1 
SUM1  -  ROl  -  TH1 
SOF  -  SIN(DIFl) 

SSM  -  SIN(SUMl) 

TOT  •  TAN (0IF1) 

TSM  -  TAN(SCMl) 

IF(SDF.LE.RMIN)  TOT  -  1.0-M 
IF (3DF. LI. RMIN)  TSM  •  1.3-M 
IF (SDF . LE . RMIN)  SSM  -  1.0-M 
IF ( SOF. LE. RMIN)  SDF  -  1.3-M 


T 


FILE:  ESCAT  FORTRAN  AL 


sii  -  (sof/ssm)««2 

RL2  -  (TDF/TSM) *«2 

FLL  -  CMPLX(2.0«(L.0-RLL) *SQRT(DL) ,0.0) 
ri2  -  CMPLX( 2 . 0* ( L . 0-RL2 ) *SQRT(DL) ,0.0) 

c 

C  INCLUDE  REFLECTION  PHASE  SHIFTS  FOR  (0,0)  RAF: 

C 

DEL  -  SQRT(ST0**2  -  K*«2)/CT0 
DELI  -  2 . 0*ATAN (DEL) 

DEL2  -  2 . 0*ATAN (N*N*OEL) 

SOL  -  SOL  -  DELL 
SO 2  -  S02  -  0EL2 
C 

C  COMPOTE  < L , 0)  AND  (0,0)  RAF  SCATTERING  AMPLITUDES  IN 

C  STANDARD  UNITS  (E.S.  VAN  DE  HULST,  WISCOMBE) 

C 

SLL  •  CMPLX(0. 0 , -L . 0*KA/2 .0) *FLL*CEXP (CMPLX( 0 . 0 , SL1) ) 

SL2  •  OfPLX(0 .0,-1. 0*KA/2 . 0) *FL2*CEX? ( CtPLX( 0 . 0 , SL2 ) ) 

60  SOL  -  CMPLX(O.O,-L.0*KA/2.O) *F0L*CSXP(OtPLX(0 . 0 , SOL) ) 

S02  -  CMPLX(0.0,-L.O«KA/2.0) •F02*CSXP(CMPLX(0. 0 ,G02) ) 

C 

C  SDK  THE  f0,0)  AND  (L.O)  RAF  SCATTERING  AMPLITUDES 

C 

SL  •  SOL  *  SLL 
S2  -  S02  -  3L2 
ANG  «  PHID 
RSL  -  REAL(SL) 

CSL  -  AIMAG(SL) 

RS2  -  REAL(S2) 

C32  -  AUtAG(32) 

IL  •  (2.0*CABS(SL)/XA)**2 
12  -  (2.0«CABS(S2)/KA) «*2 
IF(II.LE.L)  SO  TO  ill 
STOITS ( L , LLOO ) ANG , IL , 12 
L100  FORMAT (F7.4,2EL4.«) 

389  CONTINUE 

* 

C  COMPUTE  THE  OIFFRACTED  AMPLITUDES  SDL  AND  SD2 

C  WITH  CENTER  REFERENCED  PHASES 

ARG  •  KA»SIN (PHI) 

OP  ■  o.s 

IF(ARG.LT.O.O)  ARG  -  ABS(ARG) 

IF(ARG. LE.RMZN)  30  TO  90 
IF ( ABS ( PHID) .3E. 90.0)  SO  TO  99 
SO  OP  -  COS (PHI) •SJLX(ARG) 

90  SDL  •  CMPLX(DP*KA**2 ,0.0) 

99  IF(ABS(PHID) .SE.90.0)  SOL  -  CMPLX( 0 . 0 , 0 . 0) 

SOS  •  SOL 

ANG  -  PHID 

RSDL  -  REAL(SOL) 

CSDL  -  AIMAG(SDL) 

RS02  -  RSOL 
CSD2  -  CSDL 


uuooo 


FILE:  SSCAT  FORTRAN  A1 


ID  -  (2.0*CABS(S01)/KA)  *«2 
ir(JJ.lX.I)  SO  TO  389 
WRITS (1,1200) *MC, ID 
1200  FORMAT ( F7 . 4 , £1 4 . 8 ) 

389  CONTINUE 
C 

C  COMPUTE  TOTAL  AMPLITUDES  FOR  P.O.A.  PLUS  DIFFRACTION 

C 

51  -  SI  *  301 

52  «  32  -  S02 
AMO  -  PHID 
RSI  -  RSAL(Sl) 

CS1  -  AXMAS (SI) 

RS2  -  REAL(S2) 

CS2  -  AXMAC ( S  2 ) 

IT1  -  (2.0«CABS(S1)/KA) »«2 
IT2  -  ( 2 • 0*CABS (S2) /KA) **2 
IF(KX.LS.l)  SO  TO  890 
WRITS ( 1 , 13 00 ) AMO , IT1 , IT2 
1300  FORMAT (FT. 4, 2S1*. 8) 

390  CONTINUE 
C 

C  STEP  PHI  UP  BY  DPHI  AMO  REPEAT  UNTIL  PHI END 

C 

L  -  L  -  1 

PHIO  -  PHI BOM  *  L*C?HI 
IF(PHIB.LE.PHIENO)  SO  TO  50 
STOP 
END 

SUBROUTINE  FRES (WE, FS , FC) 

C 

C  COMPUTES  THE  FRESNEL  INTEGRAL  FROM  0  TO  WE 

C 

W  -  ABS (WE) 

F  -  (1.0  *  0.92S*W)/ (2.0  *  1.792*W  -  3.104*W*«2) 

S  •  1. 0/ (2 . 0  «  4 . 142*W  -  3.492«W*«2  -  3.S70«W«»3) 

A  -  1.37079«*W*«2 

C  •  COS (A) 

3  -  SIN (A) 

FC  »  0.3  *  F*S  -  S«C 
FS  -  0.3  -  F«C  -  S*S 
IF(WE.LT.O)  SO  TO  3 
FC  -  FC  *  0.3 
FS  -  FS  +  0.3 
RETURN 

3  FC  -  0.3  -  FC 
FS  -  0.3  -  FS 
RETURN 
END 

FUNCTION  BJ1X(X) 

COMPUTES  BESSEL  FUNCTION  J1(X),  AND  DIVIDES  3Y  X: 
USINS  THE  POLYNOMIAL  APPROXIMATIONS  FROM  ABRAMOWITZ 
8  STESUN 


Fill:  ESCAT 


FORTRAN  A1 


OtKEHSIOM  A(7) ,3(7) ,0(7) 

DATA  A/0. S, •0.38249919,0. 21093 97 3, -0.03984289,0. 30443319 , 

•  -3.17811-4,1.1098-3/ 

DATA  8/0. 79788438,1. 381-8, 0. 01839887, 1.71031-4, -0.00249911, 

•  0.00113933,-2.00338-4/ 

DATA  C/-2 . 33819449 ,0.12499812 , 3 . 831-3,-0 . 00837879 , 7 . 43488-4 

•  7.98241-4,-2.91888-4/ 

Z7(AJM(X)  .iT.3.0)  SO  TO  20 

Z  -  A13 ( 3 . 0/X) 

T  -  8(7) 

T  -  C(7) 

00  10  K  •  8, 1,-1 
F  -  8(X)  *  Z*F 
T  -  C(X)  •*.  Z«T 
10  CONTINUE 

8J1X  -  F*COS(XvT)/(X*SQRT(X) ) 

RETURN 

20  Z  -  (X/3 .0) **2 
8JZX  -  A(7) 

00  30  K  •  9, 1,-1 

BJ1X  -  A(K)  *  Z*8J1X 

30  continui 

RETURN 

END 


appendix  b 


SUPPLEMENTAL  Ij  EXAMPLES 

Figures  B.l  -  B.9  give  the  results  for  calculations  for  the  j  *  1  case 
(perpendicular  E-field  vector  polarization)  similar  to  those  performed  in  Chapter  2.  There 
are  several  differences  in  the  features  for  the  I{  case  which  should  be  mentioned.  First, 
the  fine  structure  variation  is  much  more  prominent  than  for  the  I2  case  as  can  be  seen  in 

Figs.  B.l  -  B.4.  This  behavior  is  also  shown  for  coated  bubble  cases.  Uncoated  bubble 
results  show  that  for  30°  <  0  <  90°  the  positions  of  the  fine-structure  intensity 
oscillations  for  j  ^  2  are  typically  shifted  from  those  for  j  -  1  by  one  half  of  a  fine- 
structure  quasi-period,  except  where  0  is  slightly  less  than  0C.  This  fine  structure  shift 
is  also  present  in  the  coated  bubble  results  although  it  appears  that  the  shift  towards  0C  in 
I]  is  slightly  greater  than  for  I2  at  comparable  coating  thickness  values. 

One  prominent  feature  seen  in  the  j  *  1  case  and  not  the  j  =  2  case  is  the 
anomalous  modulation  in  coarse  structure  for  large  values  of  h  in  the  0  -  70°  to  90° 
region  as  is  shown  in  Figs.  B.2(e)  and  B.3(e).  This  behavior  was  not  seen  in  the  j  =  2 
cases  (compare  Fig.  2.8(e)  and  2.9(e)).  The  apparent  cause  of  this  modulation  is  the 
nature  of  the  reflection  coefficient  at  the  water-coating  interface  for  the  two  cases. 
Consider  the  reflection  at  the  water-coating  interface  of  ray  (O.O)t,  shown  in  Fig.  2.7. 

From  geometric  considerations  the  local  angle  of  incidence  i  is  related  to  the  scattering 
angle  0  by 

i  =  (180°  -  0)/2.  (B.l) 

The  relative  amplitudes  of  the  waves  associated  with  this  ray  in  the  two  cases  (j  *  1  and 
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j  =  2)  are  given  by  comparison  of  the  Fresnel  (or  plane  surface)  reflection  coefficients  rj. 


sin(i  -  r)  3  tan(i  -  r) 

sin(i  +  r)  ’  tan(i  +  r) 


(B.2a,b) 


where  r  denotes  the  local  refraction  angle  and  is  given  by  r  =  sin*1  [(nw/nc)sin  i]. 
Consider  the  case  0  =  82.82°  for  which  (B.l)  gives  i  =  48.59°.  Then  for  nw  =  4/3  and 
*1.5  we  have  r-  41.81°  and  (B. 2)  gives  rj  =0.1181  and  -  -0.0008;  clearly 
tr|l »  ^l.  For  the  case  j  =  1  it  appears  that  the  amplitude  of  the  (0,0)b  ray  may  be 
sufficient  to  interfere  with  that  of  the  (0,0)  ray  in  Fig.  2.7  to  cause  the  anomalous 
modulations.  Because  this  effect  only  occurs  in  the  j  =  1  case,  the  choice  of  j  =  2  for 
bubble  sizing  and  detection  is  further  supported. 

Another  feature  which  was  much  more  prominent  in  I)  was  the  presence  of 

contrast  modulations  created  by  the  superposing  of  the  (0,0)  and  (3,1)  ray  amplitudes 
onto  the  fine-structure  interference  pattern  of  the  (0,0)  and  (2,1)  ray  amplitudes  (e.g., 

0  »  67°,  0  »  85°  in  Fig.  B.2(a)).  These  modulations  were  discussed  by  Langley  and 
Marston  for  the  case  of  an  uncoated  bubble  in  Ref.  13  of  Chapter  2.  Note  that  as  the 
thickness  of  the  coating  is  increased  other  modulations  occur  (i.e.,  0  ■  72°,  75°  in  Fig. 
B.2(d))  which  appear  to  be  unique  to  that  particular  coating  value.  This  effect  is  most 
prominent  for  smaller  bubbles  (ka  <  1000).  The  apparent  reason  for  this  behavior  is  due 
to  the  fact  that  the  approximate  angular  period  of  the  contrast  modulations  (A0)M  is 
proportional  to  X^fBj  -  B2)  where  B3  is  the  distance  between  the  (0,0)  and  (3,1)  ray 
and  B2  is  the  distance  between  the  (0,0)  and  (2,1)  ray.  The  introduction  of  a  coating 
invariably  alters  the  distance  between  B3  and  B2  resulting  in  the  new  modulations.  For 
most  other  cases  the  behavior  of  I{  was  found  to  be  similar  to  the  analogous  cases 
presented  in  Chapter  2  for  I2.  These  results  are  presented  in  Figs.  B.5  -  B.9. 


Figure  B.l.  Calculated  nonnalized  scattering  irradiances  for  ka  =  100,  nw  =  4/3,  =  1.5 

for  the  perpendicular  polarization  It  in  the  critical  angle  scattering  region.  The  solid 
curve  is  from  Mie  theory,  the  thick-dashed  curve  is  the  physical-optics  approximation,  and 
the  thin-dashed  curve  is  the  calculated  coated  sphere  results  for  (a)  0.25  pm,  (b)  0.5  pm, 
(c)  1.0  pm,  (d)  2.0  pm,  and  (e)  3.0  pm.  Note  the  greater  amount  of  fine  structure  as 
compared  to  the  I2  case  in  Figures  2.5  (a)  -  (e). 
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Figure  B.8.  Like  Figure  B.7  but  thin-dashed  curve  represents  nc  »  1.55. 
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Figure  B.9.  Calculated  normalized  irradiances  for  a  strongly  absorbing  coating 

(n,.  ■  1.5  +  ilE  -  3)  for  lea  »  1000  and  (a)  h  »  0.99  iim,  (b)  h  *  2.0  Jim,  (c)  h  »  3.14  pun, 

and  (d)  h  *  5.00  jam.  The  solid  line  is  the  non-absorbing  coat  result.  Note  that  in  case  (a) 
the  absorbing  case  (dashed)  is  virtually  the  same  as  the  non-absorbing  case. 
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APPENDIX  C 


EFFECTS  OF  THE  REFRACTIVE  INDEX  OF  WATER 
ON  SCATTERING  PATTERNS 

For  the  computations  done  in  Chapter  2  and  Appendix  B  the  value  of  nw  =  4/3 

was  used.  This  value  was  chosen  as  it  represents  roughly  the  typical  value  for  conditions 
found  both  inside  and  outside  the  lab.  Tests  were  performed  over  the  range  of  nw  = 

1.332  (pure  water  at  a  temperature  -  30°  Q  to  nw  »  1.338  (high  salinity  in  ocean  water) 
in  order  to  determine  the  range  of  validity  of  nw  -4/3.  Typical  results  are  shown  for  the 
two  extremes  in  Figs.  C.1  and  C.2.  As  can  be  seen  in  both  cases  there  is  a  slight  shift  in 
the  coarse  structure  while  fme  structure  retains  its  approximate  structure.  This  shift  was 
found  to  be  less  prominent  for  larger  bubbles.  Typically  values  of  nw  within  ±0.002 

(1.332  to  1.335)  were  found  to  have  a  negligible  effect  on  all  but  the  thickest  coating 
values  where  a  slight  amount  of  added  coarse  structure  shift  was  noted.  Under  unusual 
conditions  such  as  extreme  water  temperatures  or  high  salinity  where  variation  of  nw  is 
greater  than  0.002  the  measured  value  of  nw  given  in  the  tables  is  best  used. 

Figure  C.3  compares  the  j  -  2  imdiance  values  for  nw  -  1.332,  and  nw  -  1.338 
to  nw  *  4/3  for  a  non-coated  air  bubble  over  the  entire  180°  range.  As  can  be  seen  from 

the  figure  the  critical  angle  region  shows  good  agreement  between  all  three  values.  In  the 
majority  of  cases  the  approximation  of  nw-4/3  wu  valid  for  both  coated  and  non- 

coated  air  bubble  calculations  as  the  behavior  of  all  earlier  described  effects  was  unaltered. 
Since  most  values  of  sea  water  are  within  the  ±  0.002  range  of  nw  *  4/3  the  results 

presented  here  may  be  used  as  accurate  models  of  possible  coating  effects. 


Figure  C.I.  Water  index  of  refraction  evaluation  for  the  lower  extreme  case  of 
n  «  1.332,  (dashed  curve).  The  solid  curve  is  for  nw  ■  4/3.  The  other  parameters 


Figure  G3.  Comparison  of  inadiance  profiles  for  nw  »  4/3  (long-dashed  curve)  versus 
nw  ■  1.332  (solid  curve)  and  a*  ■  1.338  (short-dashed  curve)  for  noncoated  air  bubble  in 
water.  Note  the  excellent  agreement  between  the  three  curves  in  the  9  -  60°  to  0  •  85° 
range. 
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Following  the  completion  of  Stuart  Billette's  thesis  (which  constitutes  the  bulk  of 
the  present  report)  an  error  was  discovered  in  the  computer  program  used  to  evaluate  the 
physical-optics  approximation  of  Marston  and  Kingsbury.  This  error  (which  could  be 
traced  to  a  previous  student)  is  only  concerned  with  the  physical-optics  approximation  for 
the  scattering  from  an  uncoated  bubble. 


Furthermore, 

the  Mie  theory  results  (which  are  exact  for  an  uncoated  spherical  bubble)  were  not  affected 
by  this  error.  Indeed,  the  effect  of  the  programming  error  on  the  computed  physical-optics 
approximation  (POA)  to  the  normalized  irradiances  Ij  is  typically  only  a  few  percent  of 

the  corrected  POA  values  if  ka2  500. 

For  completeness,  the  next  few  pages  show  the  corrected  (solid  curve)  and 
erroneous  (dashed  curve)  POA  predictions  for  ka  -  100, 500, 1000,  and  2500.  The  first 
set  of  four  plots  give  I2.  The  solid  curves  of  the  plots  supersede  the  appropriate  POA 
curves  in  Fig.  2.5  and  2.8 -2.13.  The  second  set  of  four  plots  give  Ij.  The  solid  curves 
of  these  plots  supersede  the  appropriate  POA  curves  in  Fig.  B.l  -  B.4.  Except  for 
ka  » 100  the  differences  resulting  from  this  correction  would  probably  not  be  noticeable. 


The  programming  error  in  the  POA  algorithm  (FILE:  ESCAT)  is  on  page  101 .  To 
correct  the  error  change: 


CALL  FRES(W,  FC,  FS) 


to  read: 

CALL  FRES(W,  FS,  FQ. 

This  error  was  not  present  in  computer  programs  used  in  the  original  development  of  the 
POA  [Marston  and  Kingsbury,  J.  Opt.  Soc.  Am.  2 L  192-196  (1981);  Kingsbury  and 
Marston,  ibid  2L  358-361  (1981).]  Consequently  the  POA  predictions  there  are  correctly 
plotted  [see,  however,  J.  Opt  Soc.  Am.  2L  917  (1981).] 

It  is  appropriate  to  comment  on  an  additional  aspect  of  Billette's  thesis.  Recent 
tests  were  performed  which  give  additional  support  to  the  correctness  of  the  algorithm 
(FILE:  COATSPHR)  used  to  evaluate  the  exact  partial-wave  series  for  scattering  from  a 
spherical  bubble.  The  procedure  was  to  consider  the  case  of  a  shell  with  iIq  =  4/3, 
n,.  =  1.5,  and  nj  =  1  in  the  limit  a/b  -*  0  so  drat  the  radius  a  of  the  bubble  was 
negligible  in  comparison  to  both  the  wavelength  and  the  coating  radius  b.  In  this  case  the 
scattered  irradiance  went  over  to  that  of  a  uniform  sphere  (of  relative  refractive  index 
njn^)  as  predicted  from  Mie  theory. 
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